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ABSTRACT 



The reentry trajectory of a bluff zero-lift reentry body 
was computed using an HP-9845A computer. To obtain the 
external shape of an ablating reentry body as a function of 
altitude, coefficients for polynomial equations were deter- 
mined. At six discrete altitudes, the calculated nose shapes 
were identical to published information for the nose shapes 
of a reentry body with silaceous ablation material. Newtonian 
flow and an exponential atmosphere were assumed throughout. 

The computer program was modified to include ablation due 
to ice crystal erosion, and the trajectory was recalculated 
for a model ice cloud 5,000 feet thick at an altitude of 
30,000 feet. A total of four trajectories was computed for 
various values of mass loss fraction for erosion due to ice. 
The differences in range from the initial trajectory without 
ice erosion were one to fifteen feet. 
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I. INTRODUCTION 



A. BACKGROUND 

The Naval Surface Weapons Center, White Oak, Maryland, 
is. partially responsible for providing improved sensors for 
balloon-carried radiosondes which measure atmospheric proper- 
ties at various altitudes. One item of ongoing research at 
White Oak, and initially the object of this thesis, is an 
improved humidity sensor. 

An object of contention within the strategic missile 
community is whether or not water vapor, rain or ice, i.e., 
atmospheric humidity phenomena, significantly modify the 
reentry trajectory of hypersonic zero-lift reentry bodies. 

In order to determine the influence on the reentry trajectory, 
the ability to measure accurately water vapor, water droplets 
and ice crystals must exist. 

The humidity sensor currently being used on radiosondes 
is a carbon hygristor element which has poor time response 
at high altitudes because of the inherent sensing delays in 
the hygristor element. For temperatures of -40°C, which 
corresponds to altitudes above about 30,000 feet, the time 
lag of the sensing element can be of the order of fifteen 
minutes [Ref. 1]. Therefore, significant errors in humidity 
versus altitude maps are encountered. 

Table 1 is a reprint from Brousaides [Ref. 1] showing the 
time response of the currently used carbon hygristor. Various 



13 



step changes in relative humidity at several fixed temperatures 
yield the tabulated time responses to achieve 63 % or 901 of 
the step input. 



Table 1. Carbon Element Response Lags 
CReproduced from Brousaides [Ref. 1], p. 16) 



RH Level % 
(Initial- Final) 


Temp 


Change in 
63% 


RH Indication, sec 
90% 


30-50 


0°C 


0. 1 


0. 5 


51-85 




0. 3 


1.4 


85-51 




0.4 


2.2 


51-30 




0. 3 


1. 6 


31-51 


-20°C 


2. 6 


14 


50-80 




2.4 


21 


79-51 




2. 6 


48 


51-31 




3. 8 


47 


31-50 


-30°C 


15 


61 


50-72 




11 


104 


72-50 




9 


169 


50-31 




33 


163 


31-50 


-40°C 


86 


363 


50-66 




82 


27 2 


65-50 




88 


578 


50-30 




260 


7 27 



Table 1 shows that, as temperature decreases, that is, as 
altitude increases, the time response of the element becomes 
so slow as to make the readings virtually meaningless for 
any reasonable rate of upward balloon motion. In the case 
where ice is to be detected, the time lag of the sensing 
element is particularly critical. 

The object of current research on the hygristor element 
is to determine a heat transfer model of the sensing element 
so that the delay times can be calculated as a function of 
temperature. By relating delay time to temperature and 
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temperature to altitude, a more accurate map of humidity versus 
altitude can be obtained. If an accurate heat transfer model 
of the hygristor element can be obtained, data processing can 
be used to yield accurate measurements. Alternatively, an 
improved sensor can be designed once the exact cause of the 
delay can be pinpointed. 

Even if humidity can be accurately measured, determination 
of water droplet or ice crystal concentration and size distri- 
bution within a cloud is difficult. In other words, even if 
humidity can be sensed accurately, can this information be 
related directly to rain or ice? Perhaps a different type 
of sensor or several sensors might be needed to determine 
humidity, cloud composition and concentration. 

In addition, in order to design a useful sensor, criteria 
for measurement accuracy are required. The criteria for 
sensor accuracy are determined by the sensitivity of the 
reentry trajectory to changes in water vapor, rain and ice 
concentrations . 

Criteria are not available for determining the change 
in a reentry trajectory as a function of atmospheric xvater 
content. Therefore, it became necessary to start at the 
beginning in order to learn enough about the reentry process 
to make estimates of how much the reentry trajectory could 
differ with changes in humidity-related atmospheric conditions. 
This thesis estimates the change in range due to clouds with 
ice crystals. 
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Consequently, basic information concerning hypersonic 
flow about bluff bodies was studied. Equations of motion 
for reentry bodies were formulated in relation to the current 
problem. The many factors relating to reentry, such as 
ablation, convective and radiative heating and the modifi- 
cations to ablation because of chemical interaction with 
water were considered. 

B. SCOPE 

Three separate atmospheric water phenomena affect 
reentry ablation. The three phenomena are water vapor, rain 
and ice crystals. Of the three phenomena, reentry through 
ice clouds is the most destructive. Consequently, this 
thesis studies the effects of ice clouds. 
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II. THEORY 



A. REENTRY DYNAMICS 
1 . Traj ectory 

The trajectory of interest is the ballistic reentry 
of a zero-lift vehicle typical of reentry bodies used as 
warheads. A ballistic missile trajectory is typified by 
near-orbital reentry velocities and a steep initial flight 
path angle, i.e., greater than five degrees. 

All trajectories of reentry bodies pass through 
three operational phases. The three operational phases are 
defined as follows: 

a. Keplerian phase: The segment of the trajectory 

where gas dynamic forces are negligible (p = 0) . 

b. Intermediate phase: The segment of the 

trajectory where gas dynamic forces are of comparable 
magnitude to other terms (gravity) in the equations of 
motion. 

c. Gas-dynamic phase: The segment of the trajectory 

where aerodynamic forces dominate the equations of motion 

for reentry. 

Of the three phases described above, the gas dynamic 
phase is the one in which environmental factors, such as 
clouds or humidity, influence the trajectory. Therefore, 
the gas dynamic phase is the only portion of the reentry 
which will be considered. 
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The analysis of the flight of the reentry body is a 
particular application of the general theory of the dynamics 
of rigid bodies in three dimensions. The motion of the center 
of mass can be separated from the motion of the body about the 
center of mass. In this discussion, since a zero-lift body 
is hypothesized, the assumption is made that the center-line 
of the body always is parallel to the velocity vector. 

In addition, the assumption will be made that the 
reentry trajectory is planar rather than three-dimensional. 
Motion in the third dimension is, of course, possible and 
indeed certain for actual reentry bodies due to wind, 
precession, Coriolis force and so forth. However, if the 
"particle" can be described as one capable of generating drag 
only, then motion in only two dimensions is assured. Conse- 
quently, the reentry body appears as a point particle of 
varying mass subject to the various forces applied by the 
reentry environment. 

Atmospheric properties, particularly the densities 
existing at higher altitudes, are subject to considerable 
uncertainty. In order to provide a model of the density 
profile of the atmosphere, two assumptions are made: 

1. The earth and its atmosphere are spherically 
symmetrical . 

2. The atmospheric density varies exponentially 
with altitude. 

The first assumption is reasonable for this analysis 
since only a small portion of the atmosphere is traversed 
for a ballistic entry. The second assumption is based 
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on the simple kinetic theory of an isothermal gas in a 
uniform gravitational field. These assumptions lead to the 
well-known exponential model for the atmosphere: 



P = P Q e 



-By 



( 1 ) 



where: p = free stream atmospheric density (slug/ft j 
Pq = assumed mean sea level atmospheric density 
B = atmospheric density decay parameter (ft ^) 

y = altitude of vehicle above the surface of 

the earth (ft) 

B is equal to the mean molecular weight of the earth's 
atmosphere times the local acceleration due to gravity 
divided by the product of the mean atmospheric temperature 
and the universal gas constant. Figure 1 is a reprint from 
Loh [Ref. 2] and compares the exponential approximation 
with the ARDC (1956) model of the atmosphere. 

2 . Equations of Motion 

The geometry of an IRBM or ICBM trajectory is shown 
in Figure 2 [Ref. 3], Assuming no friction occurs during 
the flight, the conservation of energy gives: 



2 2 

1/2 mV - mgr = 1/2 mV^~ - mgR = constant (2) 



Equation (2) provides a unique relationship between 
velocity and altitude provided some initial condition 
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Altitude ft 1 10’ 




Figure 1. Comparison of Exponential Approximation 

with ARDC Model of Earth's Atmosphere (1956) 
(Reproduced from Loh [Ref. 2], p. 13) 




Figure 2. Diagram of ICBM Flight Path 

(Reproduced from Fuhs [Ref. 3]) 
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or boundary condition can be used to find the constant term. 
Equation (2) also can be written, if r * R, as: 




constant 



(3) 



where the constant in both equations (2) and (3) is the total 
reentry energy, i.e., kinetic plus potential energy, of the 
body . 

The conservation of angular momentum of a constant 
mass also applies and can be written as: 



Boundary conditions for equations (3) and (4) are 9(0) = 9^ 
and r(0) = R. Time can be eliminated from the equations 
using the relation from above. 



Substitution of equations (4) and (5) into equation 
(3) gives: 




(4) 



dr _ dr d9 _ dr 

<it ~ air ar d9 



RV cosA 
4 > 



(5) 




21 



(7) 



Z 



[ 



RV. cosA 
4 > 



r 



gR 

V cosA 

<P 



] 



Note that 



-R V, cosA 

dZ = dr 



( 3 ) 



Substituting equations (7) and (8) into equation (6) yields 




2gR 




(9) 



Integrate equation (9) to obtain 6 : 




( 10 ) 



If 






2 




2gR 




(ID 



then equation (10) becomes 



dZ 



( 12 ) 




arccos 




C 



To evaluate the constant C, symmetry and the boundary 

conditions r = R, 0 = ±0 and Z = Z„ can be used. Thus: 

7 o o 



z 



o 



<j> 




cosA 




ctR 

o 

cos A 



+ 




gR 



'cosA 



2 




(13) 



Substitution of the boundary conditions into equation (12) 
yields : 

Z 

cos C0 O + C) ■ COS (0 Q - C) = ~ 

This implies that C = it and therefore: 

cos (9 + ir) = - cos9 = I* (15) 

At this point, refer to equations (2) and (3), and 
note that for some value of initial velocity the total energy, 
i.e., the constant term, will be zero. This value of initial 
velocity corresponds to the initial velocity at which the 
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particle kinetic energy is equal to potential energy. Since 
the constant term is zero, as the particle moves, the constant 
term will remain zero, and kinetic energy will just balance 
the potential energy for all values of r, even as r approaches 
infinity. Thus, the value of initial velocity is the escape 
velocity. Hence 



V 



4 > 



2 




2gR 



(16) 



get 



Substitute equations (7) and (11) into equation (15) to 



COS0 



gR 

V , cosA 
JP 



RV , cosA 

<$> 



r 




- 2gR 



(17) 



Define and substitute the relations V = V/V and r' = r/R 

r 6 

and equation (16) into equation (17) to arrive finally at 



cosO 



2 - (2V 2 cos 2 A)/r ' 

2 



(1 + 4V 2 cos Z A (V 2 - 1)) 1/2 



(13) 



Equation (18) is used to calculate a family of curves 
showing range versus launch angle. In Figure 3, range is 
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found as a function of launch angle be setting r' = 1.0 in 
equation (18) to find 0 Q for various values of V between 
0.5 and 0.7 and realizing that 

Range = 2R0 Q 

7 

where R = 2.09 x 10 ft for the radius of the earth. 




UinCH i BGL2 



Figure 3. Range of a Ballistic Missile as a Function of 
Launch Angle for Various Values of V = V /V 
(Reproduced from Fuhs [Ref. 3]) ® e 



( 19 ) 
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A value of can be deduced from Figure 3 by assuming 
that the missile will be launched with the proper velocity 
and launch angle to achieve the optimum range, i.e., the peaks 
of the curves. 

For example. Table 2 shows the various characteristics 
of typical ballistic missiles in the inventory today. The C-4 
Trident has a range of between 4000 and 4800 nm. Referring to 
Figure 4, an optimum range of 4100 nautical miles requires a 

launch angle of about 27° and V of about 22006 ft/sec. 

<P 



Table 2. Launch Vehicles and Parameters 
(Reproduced from Fuhs [Ref. 3]) 



Missile 


Ola 

Icfiqln 


Weight 


Payload 


P.inge 

mo 


Warhead 

Yield 


10C 


Re 1 1 red 


CEP 

feet 


Stages 


A-1 Polaris 


54-2fl* 


30,000 




1200 




Nov 60 


Oct 65 




2 


A- 2 Polaris 


54"3I ' 






1500 










2 


A-3 Polaris 


54**31 ' 






2500 










2 


C-3 Poseidon 


72'* 34 * 




14RV 


2300 








1800 




C-4 Trident 


72" 34* 


70,000 


8MK4 


4000 
(4800 ) 


100KT 


1979? 




1500 


3 


0-4 Trident 


S?“45. 8* 






6 700 












0-5 Trident 


132-45 .8* 


126,000 


14P.V 

_JMK17A 




350JCT 


1005 

(1937) 








MX 


92 M 


172,000 


KMK12A 


6000 




1983 




(100* dne to 
qui dance ) 
?30' r^o) 


Atlas 














1965 






Titan 






MK6 














MM I 




















rein 




















Mill 






3MK12A 


B000 


350KT 






600 


3 
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Returning to equation (2) , the constant for the C-4 



Trident can be found by substituting and R. 



1 y 2 _ p = constant , 
1 <p S ~ m 



( 20 ) 



K = j (22006) 2 - (-32.2) (2.09 x 10 7 ) = 9.15 x 10 8 ft 2 /sec 2 



Since K now is known, equation (2) can be used to find the 
velocity at any point in the trajectory given r by rewriting 
equation (2) as 



V = [2(K + gr) ] 1/2 



( 21 ) 



If the portion of the trajectory near reentry is 
examined (that is, where r is not much different from R) , 
the earth can be approximated as being flat since y/ft << 1 
where y is altitude. For this case, the simple two- dimensional 
equations of motion in Cartesian coordinates apply, that is 

IF = mx (22) 

X 

and 

EF y = my (23) 
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where x is distance measured horizontally, and y is distance 
measured vertically. 

Since in the model the only force in the x-direction 
is drag, and the forces in the y-direction are drag and 
gravity, equations (22) and (23) can be written as 



x = 




(24) 



y 




m 



(25) 



where D and D represent the drag in the x and y directions, 
x y 

respectively. In finite difference form, equations (24) and 
(25) are 



X 



3 + 2 



2X . + 

]+l 



At' 



X. D 
3 = x 

m 



(26) 



Y . - 2Y . + Y . 

3+2 3+1 3 = . 

I? 



D 



g - 



y 



m 



(27) 



The values of X^ and Y^ and V(o) are obtained from initial 
conditions at the point of interest at the beginning of 
reentry using equations (18), (19), and (21). The values of 
X 2 and Y 7 are obtained from 
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X, + U At 
1 o 



( 28 ) 



X 



2 



Y = Y + V At 
2 1 o 



( 29 ) 



where 




( 30 ) 



B. HYPERSONIC REENTRY 

The character of the gas dynamic phase of flight at 
near-orbital velocities is far removed from that of flight at 
sonic or subsonic velocities. The vehicle has tremendous 
energy and angular momentum, most of which must be released 
to the atmosphere. The two processes available to transfer 
the energy are manifested as heating of the vehicle and 
deceleration . 

No precise definition may be given for the velocity 
at which supersonic flow becomes hypersonic flow because 
the effects vary with the vehicle geometry. In addition, the 
onset of these effects occurs gradually rather than abruptly. 
Knowledge of Mach number at which transition from supersonic 
to hypersonic flow occurs is important. If the flow is hyper- 
sonic, certain simplifying assumptions can be made concerning 
the shock shape and position relative to the body, which 
relate directly to drag- and pressure coefficient. These 
assumptions are discussed below. 
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During hypersonic flight, there will be regions of the 
flow fields about a typically blunt-nosed body where the 
flow decelerates and in which subsonic, transonic and 
supersonic flows all exist simultaneously. For the purposes 
of this discussion, a Mach number greater than or equal to 
four will be considered hypersonic. 

Any bluntness of the nose will manifest itself as a 
detached shock. A pointed nose will melt or ablate due to 
the intense heat generated at the point and therefore will 
become blunt. Consequently, the blunted nose generates the 
detached shock. The most striking effect of hypersonic flight 
is the formation of a shock ahead of virtually all bodies. 

The behavior of the shock for the flow past a blunt body of 
revolution as the Mach number increases may be illustrated 
as in Figure 4, taken from Cox and Crabtree [Ref. 4], for 
flow of a perfect gas about a sphere. 

As can be seen, the shock stands off from the body at all 
Mach numbers, and the shock-detachment distance decreases 
as the Mach number increases. The limiting shock (M -> °°) lies 
close to the body and is very nearly concentric over most of 
the surface of the sphere. Just behind the shock wave where 
the shock is nearly perpendicular, the flow is brought almost 
to rest and is subsonic. As the flow passes around the body, 
the flow accelerates, passes through sonic speed and finally 
accelerates to supersonic velocities as is illustrated in 
Figure 5, again taken from Cox and Crabtree [Ref. 4]. 
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Figure 4. Shock Shapes for Flow Past a Sphere at Different 
Mach Numbers (y = 1.4). (Reproduced from Cox 
and Crabtree [Ref. 4], p. 9) 

At the high temperatures experienced behind the shock in 
hypersonic flow, the assumption of a perfect gas cannot be 
made as the energies of the gas associated with temperature 
are comparable to the energies of various molecular and atomic 
processes such as ionization, dissociation and vibration. 
Therefore, characteristics of the gas such as the commonly 
accepted value of the ratio of specific heats for air, i.e., 

Y = 1.4, must be viewed with skepticism and taken into account 
when calculating the flow field. 

However, the Mach number can be associated with the ratio 
of the kinetic energy and the random thermal energy by the 
relation 
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(31) 



M 



2 



(velocity) V 2 
(speed o£ sound) /2 



* v!zz . « 

yRT/2 



kinetic energy 
thermal energy 




Figure 5. Subsonic and Supersonic Regions for Flow Past a 
Blunt Body. (Reproduced from Cox and Crabtree 
[Ref. 4] , p. 65) 

For large Mach numbers, regardless of the value of y, the 
kinetic energy is much larger than the thermal energy. This 
implies that the random motion associated with collisions 
between atoms and molecules is of relatively little importance 
compared to the kinetic energy of the particles and therefore • 
the gas may be thought of as a stream of non-interacting 
particles. The stream of non-interacting particles was first 



32 



conceived by Newton for a "rare" gas composed of particles with 
large distances between them. Newton showed that if the 
particles do not interact with one another and at large 
velocities, then the resistance of a body into which the 
particles collide will vary as the square of the velocity. 

In addition, he calculated, by assuming that there was no 
transfer of tangential momentum and that the normal momentum 
of the particles is destroyed at the surface, that the force 
on a plate at incidence to the flow varies as the square of 
the angle of incidence. More recently these assumptions have 
been used to show that the pressure on the plate is given by 



Also, for flow over wedges at high Mach number, the 
assumptions were used to show that the relation between the 
shock angle (u) and the body angle (ip) is 

U _ Y + 1 n 

V 2 

Therefore as y + 1, the body and the shock coincide and 
the density ratio across the shock becomes infinite. Thus 
there is a layer of high density gas lying near the surface. 



AP = pV“sin^<p 



( 32 ) 



or 



C 



P 




( 33 ) 
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This is the "Newtonian Shock Layer" or entropy layer since 
very large changes in entropy also occur. 

Lees [Ref. 5] adjusted the Newtonian Formula [equation (33)] 
empirically so that the value at the stagnation point on the 
nose is exact. Lees' result, which is called the "Modified 
Newtonian Formula," is given by 



C 



P 




. 2 
sm ip 



(35) 



where C is the stagnation point pressure coefficient. The 
*o 

value of C is given by the Rayleigh formula developed from 
P o 

the integral Bernoulli equations and the shock equations and 
is 




2 

yM 2 



Y + 1 
2yM Z - y 



+ 





1 



(36) 



As in the Newtonian formula, as M + » and y -*■ 1 , the value 

of C approaches two. 

Po 

There are various other modifications to the Newtonian 
Formula including the Busemann Formula, the Modified 
Newtonian -Bus emann Formula and others but, as is shown in 
Figures 6 and 7, a comparison of the surface pressure 
coefficients for the flows over a sphere and a cylinder at 
M -*■ «> and Y = 1.4 confirms that the Modified Newtonian 



34 




Figure 6. Surface Pressure Distributions for Flow Past a 
Sphere (M = 00 , y = 1.4). 

(Reproduced from Cox and Crabtree [Ref. 4], 
p. 73). 
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Figure 7. Surface Pressure Distributions for Flow Past a 
Circular Cylinder (M = y = 1.4). 

CReproduced from Cox and Crabtree [Ref. 4] , 

P- 72). 
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Formula is in best agreement with the actual numerical 
calculations . 

1 . Aerodynamic Heating and Ablation 

During hypersonic reentry, excessive local heating 
is a serious problem. For the ballistic vehicle, where the 
lift-to-drag ratio is zero, very high heat fluxes are 
encountered but for a relatively short duration, i.e., five 
to 100 seconds. In practice, the heat and mass transfer 
phenomena associated with this reentry mode are nonlinear 
functions of time due to the nonlinear environmental parameters, 
chemical reactivity, mechanical and acoustical vibration, 
gas temperature, mode of energy transport, total heat load, 
peak heating rate, shape of the heat pulse, particle impinge- 
ment, material properties, etc. Of primary concern is the 
means of accommodating the energy converted to heat. 

The thermal protection of the reentry body can be 
accomplished in three ways: (a) absorption, (b) radiation, 

and (c) reduction of aerodynamic heating prior to reaching 
the vehicle. The added requirement to minimize body weight 
in order to achieve greater missile range has often led to 
the choice of materials designed to change phase when heated 
and to vaporize at the surface exposed to the hypersonic 
environment. The energy absorption achieved through such 
processes has proved to be highly effective. 

The energy transfer mechanisms between the boundary 
layer gas and the ablating surface are shown in Figure 8, 
reprinted from Hurwicz, Kratsch, and Rogan [Ref. 6], The 
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Figure 8. Heat Flux Balance at the Surface. 
(Reproduced from Hurwicz, Kratsch, and Rogan [Ref. 6], p. 9) 

energy balance at the surface established the boundary con- 
dition for solution of the boundary layer equations. The net 
heat flux through the surface to the interior of the ablator 
(q s ) equals the summation of the aerodynamic heat flux to 
the surface from the boundary layer (-q ) , the energy trans- 
fer away from the surface by normal gas mass flow (m^h^) , the 
absorption of incident radiation ( a w 4 rac j) , the reradiation of 
energy from the surface ( + e w crT w 4 ) , the energy transfer by 
mechanical erosion (ablation) (m r h s ) , the energy transfer to 

the surface by condensed material (m h ) , and the energy 

c s c s 

transfer to the surface by internal gas products (m CT h ) . 

S s S s 

Note that the aerodynamic heat flux (q ) is the combination 

W 
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of convection to the surface as well as chemical energy 
transfer due to surface chemical reactions from mass transfer. 
The preceding factors can be expressed as an energy equation 



q s = 



x w 



- in h + 
w w 



m 



c s c s 



+m h - e crT + a q , - m h (37) 
g c g c w w w n rad r s v J 



Equation (37) is applicable for both laminar and 
turbulent flow and couples the energy balance to the boundary 
layer solution which is required to obtain species concentra- 
tions and the heat and mass transfer analysis. 

Equation (37), although general, does assume thermo- 
chemical equilibrium at the surface and in the boundary 
layer. Such equilibrium may not exist in reality since the 
mass ablation rate probably is reaction rate controlled by 
the ablation material, its char products, and the boundary 
layer species. In this case, gaseous ablation products can 
be determined in principle through the simultaneous solution 
of: (a) all surface thermochemical equilibrium reactions 

expressed in partial pressures and equilibrium constants, 

(b) equations of conservation of elemental species at the 
surface, and (c) Dalton's law of partial pressures. 

Even though conceptually the surface thermochemical 
equilibrium reactions can be expressed explicitly, the actual 
products of the reactions are difficult to determine because 
of the elevated temperatures and the number of interacting 
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elements. If, for example, carbon is used as an ablator, 
the products normally considered are C, C~ , , C^, C^, CO, 

C0 ? , CN, 0 7 , . As can be seen from Figure 9, the concen- 

trations of the gases change as temperature (or velocity) 
changes. In addition, if water vapor is included, as is of 
interest in this discussion, another element, hydrogen, must 
be included as a possible reactant. Thus, entire organic 
chains of molecules may be formed. 




Figure 9. Mole Fraction Concentration behind Normal Shock 

(Reproduced from Cox and Crabtree [Ref. 4], p. 10) 

In view of the complexity of ablation with non- 
equilibrium chemical compositions, information concerning 
nose shape was sought in lieu of direct calculation. Informa- 
tion was found relating nose shape to specific altitudes. The 
nose tip chosen was composed of a silaceous material. 
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Ablation occurred during reentry through the atmosphere 
including mechanical erosion. Further, turbulent floxv was 
assumed to exist from the start of reentry and continue 
throughout the flight. Figure 10 shows the family of nose 
shapes selected for this analysis. The material loss due to 
ablation is excessive compared to modern carbon nose tips; 
however, the effects of nose shape change on the trajectory 
were to be examined. Therefore, the larger the ablation 
rate, the more the trajectory could be expected to change 
from the non-ablating entry. 




Figure 10. Nose Shape of Reentry Body at Various Altitudes. 

(Reproduced from Hurwitz, Kratsch, and Rogan 
[Ref. 6] , p. 16) . 



Recent experiments have shown that interesting after- 
body ablation patterns are formed during reentry. 



41 



Cross-hatching or other similar spiral ablation patterns 
occur, and related roll torques are induced producing erratic 
perturbations in roll behavior. Such phenomena were not 
considered in this analysis. 

2 . Drag 



coefficient, ip is the local slope of the body, and A is an 
element of area. Given the nose shapes in Figure 10 and 
assuming that the shape of the nose can be determined for 
every altitude, the drag at any altitude can be computed 
using equations (35) and (36) . For an axisymmetric body 
where r = f(x), the body can be divided into infinitesimal 
segments of X; and the drag can be found on each elemental 
circular element. Figure 11 shows a typical element in the 
x-r plane. Figure 12 shows the same element as seen along 
the longitudinal axis and where R is the base radius of the 
reentry body. Therefore, the incremental area is 



The drag of a reentry body experiencing Newtonian 
flow can be expressed as : 




A 



(38) 





( 39 ) 
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Figure 11. Cross-Section of Elemental Area along the Nose. 




Figure 12. Front View of Elemental Area along the Nose. 
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Hence, equation (38) can be written for the elemental drag 
and becomes 



d (Drag) 




r + dr 

Cp 2nrdr 



( 40 ) 



The total drag on the body is then the sum of the 
elemental drags over the length. 



x=L 

Drag = \ qCpir ( r 2 2 - i^ 2 ) (41) 

x=o 



Since the Newtonian flow model applies over the 
length of the nose, the "afterbody'' can be included in the 
nose shape equations if properly written; hence equations (38) 
to (40) will define the drag in all regions along the body. 

For projectiles which fly at low supersonic velocities, 
the base drag is very important when considering the total 
drag on the body. However, at hypersonic velocities, the 
pressure increase across the shock at the nose is so large, 
of the order of one hundred atmospheres , that the low pressure 
at the base of the reentry body is comparatively insignificant. 
In this analysis, the base pressure drag was neglected. 
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C. HUMIDITY EFFECTS 



As noted previously, the effects of water vapor, water 
droplets and ice crystals in the atmosphere lead to extremely 
complex modifications to all the equations previously cited. 
Considering water vapor only, density of the air is changed. 

As a result, the modified density influences reentry body 
velocity, ablation rate, etc. Also, as mentioned before, more 
complex chemical equations are required to determine ablation 
products and, hence, heat loss and absorption. The total 
effect on drag, for instance, would be the sum of the effects 
on each parameter affecting drag; for example: 



Drag _ 
Drag 



A density 
density 



AC 

_J2. 



Aq 

q 



A area 
area 



( 42 ) 



If rain, i.e., suspended water droplets, also is included 
in the analysis, mechanical erosion must be considered. 
Additionally, when the droplets cross the shock wave in front 
of the vehicle, scattering and annihilation will occur; and 
only some fraction of the droplet mass actually will contact 
the surface of the body. The interaction of water vapor and 
rain with the reentry body is not considered in this thesis. 

The interaction of ice in high altitude clouds would seem 
to be a worse case of humidity effect interaction with a reentry 
vehicle since the assumption can be made that most of the 
typical ice particles will survive after passing through the 
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shock front. Hence, mechanical erosion should be increased 
greatly. Experiments have shown that the mass ablated can 
be many times the impact mass of the typical ice particle. 

Equations also are simplified since the time in the 
boundary layer is so short, of the order of microseconds, 
that the ice will not undergo a change in phase before 
colliding xvith the body. Therefore, changes in ice crystal 
properties and chemical interactions were neglected. 

Since information was not available concerning the 
actual change in nose shape due to collisions with ice 
particles, the assumption was made that ablation shape was 
not changed but merely was increased in rate of mass loss. 
That is, the identical equations describing the nose shape 
were used in both the baseline trajectory model and the ice 
cloud trajectory model. However, as mass loss increased 
due to impact, the nose shape was changed according to those 
equations until the proper mass loss was achieved. This 
assumption is no doubt incorrect in that it would be expected 
that the mass loss would be much greater near the stagnation 
point and much less along the perimeter of the body. 
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III. SIMULATION 



A. BASELINE TRAJECTORY 

In order to evaluate the consequences of reentry through 
ice clouds, a reference trajectory was needed to compare such 

i 

parameters as range, impact velocity, changes in drag coeffi- 
cient, etc. The equations of motion were presented previously, 
but initial conditions and boundary conditions are required 
to solve for a specific trajectory. 

1 . Selection of a Typical Reentry Vehicle 

The model selected as the reentry vehicle conforms 
to the nose shapes shown in Figure 10. In addition, realistic 
values of ballistic coefficient and drag coefficient were 
required. Modern ballistic reentry vehicles have ballistic 
coefficients (weight divided by the quantity drag coefficient 
times base area) in the range of 1500-2000 pounds per square 
foot and drag coefficients slightly greater than 0.1. 

The weight of the reentry vehicle was not known, 
but a specific gravity of 1.4 was selected as typical. 

During the course of the initial iterations of the computer 
model, the body length was changed to add enough mass so that 
the ballistic coefficient was near the desired value of 2000. 
The body shape was achieved by extending the after-section of 
the initial nose shape shown in Figure 10. Figure 13 shows 
the final model shape prior to reentry. 
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Figure 13. Model Reentry Vehicle Shape 



2. Initial Conditions 



Referring to Figure 4, the center trajectory was 

chosen as a representative trajectory of a typical Navy 

reentry vehicle which might be launched by Polaris or Trident. 

This curve corresponds to a maximum range of about 4100 

nautical miles if launched at an angle of 27° with an initial 

velocity of 22006 feet per second. Substituting the values 

for launch angle and initial velocity, equation (20) yields 
8 2 2 

K = 9.15 x 10 ft /sec as shown previously. 

Since deceleration and ablation have not yet begun 
at 250,000 feet, this altitude was chosen as the initial 
altitude to begin reentry. For an altitude of 250,000 feet, 
the value of r, the distance to the center of the earth, is 

7 

2.115 x 10 feet; the value for r was substituted into 
equation (21) to yield 

= [2 (K + gr)] 1/2 = 2 [ 9.15 x 10 8 - 32.2(2.115 x 10 7 ) ] 1/2 

= 21,632 ft/sec (43) 

At an altitude of 250,000 feet atmospheric density is too low 
to cause drag. Also, at that altitude r is not much different 
from R, the radius of the earth. Therefore, the coordinate 
system can be shifted from r,9 to x,y. Since x forces are 
zero, the x component of velocity will remain constant and 
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will be equal to the launch velocity times the cosine of the 
launch angle 



V = V cosA = 22006 cos 27° = 19,607 ft/sec (44) 

X (J) 

At 250,000 feet, the x-component of velocity must be 
the same value; and since the value of V is known, the 
reentry angle A can be determined from 



A 



cos 



-1 




cos 



-1 19607 _ 
21632 



25.04° 



(45) 



Checking these values with the curve selected in 
Figure 10 and using equation (18) , 



cos 0 Q 



1 - 2 ( ♦ 6 ) 2 cos 2 (27°)/l 

[1 + 4(.06) 2 cos Z (27°) (,6“ - 1 )] 1/2 



0.827 



Therefore, the initial angular distance from apogee is 
0 Q = 34.21°. Thus total range using equation (19) is 



Range = 2R9 0 



„ 2 (2.09 x IQ 7 ) (34.21)* 
6000 ^31) 



4160 nm 



which corresponds very well to the expected value of 4100 nm. 
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The range the reentry vehicle is expected to fly can 
be found from the reentry location to impact. To find the 
range for a no-drag trajectory recall that the arc length 
traveled is 



9’ = 0 - 9 

o 



( 46 ) 



Again, using equation (18) and substituting values for 
V, A, and r' at the entry point 



_ 1 - 2(.6) 2 cos 2 (25.04°) / (1.01) 

cos 9 - * j ^ 7 T7T 

[1 + 4(.6) cos (25.04°)(.6 - 1)] ' 

and 9 = cos " 1 ( . 843) = 32.54°. 

The angular arc of flight using equation (46) is 



0.843 



0. = 9 q _ 9 = 34.21° - 32.54° = 1.67° 
and using equation (19) 



Range = R9 ' 



(2.09 x 10 7 ) (1 . 67)77 
6000 180 



101.5 nm 
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Since the sequence of equations above assumed no drag, the 
range should be the maximum expected range of the model 
reentry vehicle . 



Table 3 



Summary of Initial Conditions at 
an Altitude of 250,000 Feet 



— 1 


Initial 


Velocity 


21,632 ft/ sec 1 


Initial 


Entry Angle 


25° 


Initial 


x velocity 


19,607 ft/sec 


Initial 


y velocity 


9,138 ft/sec 


Maximum 


range 


101.5 nm 



3 . Nose Shape 

The nose shape is the primary factor in determining 
total vehicle drag. Since the nose shape was given only for 
discrete values of altitude as in Figure 10, the nose shape 
was constructed artificially as a continuous function of 
altitude. This was done by using the HP-9845A computer along 
with the HP-9874A Digitizer. Each of the six nose shapes was 
divided into two segments ; one segment was from the nose tip 
to 2.4 inches aft of the nose tip, and one segment was from 
2.4 inches aft of the nose tip to about 20 inches aft of the 
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nose tip. Another function was required to describe the 
recession of the nose tip. 



Each nose shape was traced with the digitizer 



inputting data points into the computer. A curve fitting 
routine was used to construct curves of nose radius as a 
function of length along the centerline of the reentry body 
using the data from the digitizer. The best fit curve for 
the segments from the nose tip to 2.4 inches aft of the nose 
tip was a fifth-order polynomial with constant coefficients. 
The segment from 2.4 inchds aft of the nose tip was described 
by a straight line. Appendix A contains graphs of each of 
the six nose shapes (two segments each) comparing the data 
input with the plot of the equations which describe the data. 
As can be seen, the graphs very closely approximate the 
actual nose shapes. 



The nose shapes thus were defined mathematically as 



r 



i=o 




(47) 



for 0 < x < 2.4 inches and as 



1 

r h £ b^x 



i=o 



l 



(48) 
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for x > 2.4 inches where r is the local radius of the 
reentry body and x is the distance along the centerline. 



and (48) , plots of each coefficient versus altitude were 
constructed. For example, all the fifth-order coefficients 
of equation (47) were plotted at the altitudes applicable 
to the particular nose shapes. The coefficients in equation 
(47) thus were defined by 



where h is altitude in feet. The coefficients in equation 
(48) were defined by 



Division of the coefficient curves into two portions was 
necessary; equations (49) and (50), sixth-order polynomials, 
applied below 90,000 feet. Equations (51) and (52) 



Using the coefficients calculated in equations (47) 




(49) 




(50) 




( 51 ) 
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B. 

l 



( 52 ) 



1 

E 

j=o 



F- -h 
iJ 



j 



were used above 90,000 feet. 

Equations (47) through (52) can be combined to 
yield the nose shape, the radius, as a function of altitude 
in the various regions of applicability. Table 4 shows the 
combined equations (53) through (56) and the regions of 
applicability . 



Table 4 

Summary of Nose Shape Equations 



Altitude 


0 < x < 2.4 (inches) 


x > 2.4 (inches) 


Equal to 
or Below 
90,000 
feet 


r = Z Z C,,h J x 1 (53) 

i-0 j-0 ^ 

! 


1 

16 j i 

r = Z E D- -h J x (54) 

i=0 j =0 ^ 

, 


Above 

90,000 

feet 


51 i i 

r = E E E..h J x x (55) 

i=0 j =0 ^ 


11 i i 

r = E E F, .h-V (56) 

i-0 j =0 1J 



Appendix B compares graphs of the polynomials 
describing the coefficients versus altitude with corre- 
sponding graphs of the digitizer data input used. Appendix C 
contains similar graphs of the nose tip recession versus 
altitude. All graphs show that the equations adequately 
represent the described data. 



55 



The result of the above analysis was to allow the 
insertion of any altitude into equations (53) through (56) , 
yielding the corresponding nose radius as a function of the 
centerline distance from the nose tip. 

4 . Flow Chart 

Using the initial entry conditions and nose shape 
versus altitude information previously described along with 
the equations of motion, a program to calculate trajectory 
was possible. Figure 14 describes the process. 

5 . Validity 

The program initially was modified to print the nose 
shapes every second along the trajectory so that a comparison 
between the desired nose shapes described in Figure 10 and 
the nose shapes described by the equations derived in 
Section IIIA3 could be made. Figure 15 shows a family of 
curves plotted from the program output for nose shape at 
altitudes closely corresponding to the altitudes for the family 
of curves in Figure 10. Although the nose shapes in Figure 15 
are somewhat more square in appearance than those in Figure 10, 
the proper general shapes were clearly achieved with the nose 
shape as a function of altitude equations. 

A goal of the model was to achieve a drag coefficient 
of about 0.1 and a ballistic coefficit of about 2000. In the 
simulation, the drag coefficient ranged from 0.100 to 0.127 
during the flight with the initial value of 0.105. The 
ballistic coefficient ranged from 1392 to 1878 during the 
flight with the initial value of 1798. The range should be 
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Figure 14. Flow Chart of Baseline Trajectory Program 



57 



(Note: Nose shapes progressing from left to right correspond 

to altitudes progressing from top to bottom.) 
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figure 15. Nose Shapes at Various Altitudes 



less than the vacuum trajectory range of 101.5 nm; the range 
with drag was 84.103 nm. 

For ballistic entries at escape velocity and initial 
angles of inclination greater than about 5°, the variation of 
Velocity with altitude has been computed by Loh [Ref. 7] for 
various values of ballistic coefficients. Figure 16 shows a 
plot of the various deceleration curves. A similar curve can 
be generated from the simulated reentry. The generated curve 
should be the same shape as Figure 16 but since the reentry 
velocity is only about 80 percent of escape velocity, all 
values of dimensionless velocity will be lower. 

In addition, the simulated reentry body undergoes 
ablation and hence changes the value of the ballistic 
coefficient. However, if a nominal ballistic coefficient 
of 1500 is chosen, the value of 



C D A 1 a. 1 

WsinQ " 3sin9 1500 sin 25 u 



0.0016 



(57) 



can be us,ed, where 9 corresponds to reentry angle. Quantita- 
tively, the simulated reentry curve should lie between the 
0.001 and 0.01 curves and to the left of the curves in 
Figure 16. 

Figure 17 is a graph of the Mach number versus 
altitude for the model. Maximum deceleration occurs at 
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400 



Figure 16. 




Velocity Variation during Direct Penetration of 
the Earth’s Atmosphere. 

([Reproduced from Loh [Ref. 7], p. 110) 
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Figure 17. Altitude versus Mach Number 



about 75,000 feet, and the actual impact Mach number is 4.42. 
Based on the above parameters and the similarity in shape 
between curves in Figures 16 and 17, the computer model was 
accepted as a valid representation of the actual reentry 
trajectory . 

B. TRAJECTORY THROUGH AN ICE CLOUD 
1 . Model Cloud 

In order to evaluate realistically the effects of 
ice clouds on the trajectory of a reentry body, a repre- 
sentative ice cloud was needed. Such information as altitude, 
thickness, particle size and concentration was needed. 

Although information such as this is readily available for 
relatively low altitude rain clouds , useful information 
regarding high altitude ice clouds is very rare. The only 
information found is contained in Table 5, reproduced from 
Heymsfield [Ref. 8], in which aircraft measurements of five 
different cirrus uncinus clouds were made. 



Table 5 

Summary of Cloud Measurements 
(Reproduced from Heymsfield [Ref. 8], p. 800) 



Calculated Mean.) 

precipita- crys- I 

tion rate* tal 

(mm hr -1 ) (mm) 

length 



14 


*9 July 1070 


Ci unc 


6.5-8. 8 


6. 5-8.0 


-15 to —JO 


0.30 


0.24 





0.019 


17 


9 


1.0 


0.6 


0.9 


CM) 


10 Sent 1071 


Ci unc 


8.5-0. 1 


8.5-9. 1 


-.16 to -4.1 


0.30 


0.24 


— 


0.025 


36 


14 


1.6 


0.7 


0.75 


C(2) 


10 Sept 1071’ 


Ci unc 


8.0-0. .15 


8. 0-0. .15 


40 to -44.5 


0.25 


0.15 


— 


0.03 


14 


7 


0.97 


0.5 • 


0.75 


G(I) 


6 Apr 1072 


Ci unc 


0.6-1 1.0 


0.6-1 1.0 


-41.7 to -57.6 
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The table shows that the clouds ranged in altitude 



from 4.45 to 11 kilometers (14,600 to 36,100 feet) and 
varied in thickness from 0.3 to 3.45 kilometers (980 to 
11,300 feet). Based upon the above figures and the advice 
of Dr. K. L. Davidson, Meteorology Department, Naval 
Postgraduate School, the representative ice cloud was 
selected to be at an altitude of 30,000 feet with a thickness 
of 5,000 feet. Ice particle concentration at that altitude 
nominally was chosen as 0.24 gram per cubic meter which 
corresponds to 4.65 x 10 ^ slug per cubic foot. 

2 . Mass Defect 

The mass of the reentry body ablated by collision 
with the ice particles can be expressed by 



Mass Defect 



Base area x Cloud thickness x Mass Loss Factor 

cQS A 



(58) 



where the base area of the reentry body and the angle of the 
trajectory with respect to the horizontal, A, are computed 
at the proper altitude by the computer program. The Mass 
Loss Factor represents the ratio of the mass ablated to the 
mass of the ice particles which collide with the body. Some 
information for Mass Loss Factor exists for multiple impacts 
without a shock wave ahead of the reentry body, but none could 
be found which included the shock wave. Therefore, several 
nominal values of Mass Loss Factor were selected, and the 



matrix of 10, 20, 50, and 40 was computed. A trajectory 
was calculated for each. 

One assumption inherent in equation (53) is that all 
of the ice particles contained in a tube representing the 
base area of the reentry body passed at the proper angle 
through the cloud thickness do impact on the body. Such an 
assumption was made since this analysis was based on a worst 
case scenario. If only a portion of the ice impacts on the 
body, which is the more likely case, less ablation will occur 
and, therefore, less of an effect will result. 

3 . Flow Chart 

Using the cloud information contained in Table 5 
and equation (58) , the flow chart in Figure 14 was modified 
to include ablation due to ice erosion. Figure 18 describes 
the process . 
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Figure 18. Flow Chart of Trajectory Program with Ice Cloud 
Ablation 



IV. RESULTS , CONCLUSIONS AND RECOMMENDATIONS 



A. RESULTS 

Table 6 compares the difference in range of the various 
trajectories modeled as functions of Mass Loss Factor, which 
is the ratio of mass loss of ablation material due to erosion 
to the mass of ice particles impacting. Also shown are the 
values of mass loss due to ice erosion and final Mach number. 
Figure 19 shows the same results graphically. 

Table 6 

Comparison .of Trajectories with Ice Erosion 
to Baseline Trajectory 



j 



Mass Loss Factor 


Final* Mach Number 


Mass Loss Due to Ice (elu*a) 


Difference in 3ange (feet) s 


0 


4.417 


0 


0 


10 


4.418 


0.0947 


1.54 


20 


4.420 


0.1895 


2.85 


30 


4.425 


0.2842 


7.99 


40 


4.429 


0.3789 


15.05 



B. CONCLUSIONS 

The most interesting result of the analysis is that the 
reentry vehicle with ice erosion traveled farther in all 
cases than did the vehicle without ice erosion. The reason 
is that, despite the significant mass loss due to erosion 
which tends to decrease the ballistic coefficient, the final 
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MRSS LOSS FRCTOR 



nose shape, as seen in Figures 10 and 15, reduces the drag 
coefficient and tends to increase the ballistic coefficient. 

The result is that the two parameters, mass loss and drag 
coefficient, are of opposite sign and tend to cancel. 

Circular Error Probability, CEP, is a constant multi- 
plied by the variance in range, and for a Gaussian distribution 

a 2 = E a. 2 (59) 

i 1 

where o^ denotes the variance in range caused by the various 
factors that influence changes in range. The assumption can 
be made that another ck , that is, the variance due to ice 
erosion, can be added to equation (59) and further that the 
variance due to ice erosion is of the order of the change in 
range due to ice erosion. If it also is assumed that the CEP 
is of the order of the variance, then equation (59) can be 
written 



CEP 2 = E CEP i 2 = CEP^ + A Range 2 (60) 

where CEP^ is the original CEP of the body due to all 
factors excluding ice erosion and A Range is the difference 
in range due to ice erosion. For a C-4 Trident with an 
initial CEP of 1500 feet (from Table 2) and a Mass Loss 
Factor of 40 for a A Range of 15.03 feet, equation (60) 
yields a new CEP of 1500.08 feet. 
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The assumptions in this thesis were made on the basis of 
the best information available in the literature. The analysis 
of the problem based on the assumptions indicates that the 
difference in CEP due to ice erosion is very small indeed. 
Considering ice erosion as the worst case for atmospheric 
water phenomena, the CEP change due to water droplet erosion 
or due to differences in humidity also can be considered to 
be very small. 

C. RECOMMENDATIONS FOR FUTURE WORK 

1. The nose shape after impact with ice particles is 
obviously very important. Experiments which measure the 
effects of multiple collisions of ice on nose shape are 
required. Such experiments must include the effects of a 
bow shock around the nose equivalent to a Mach number of 
about 10. Conceptually, more ice should impact near the 
centerline of the vehicle and less toward the edges, tending 
to make the body more blunt and therefore tending to increase 
the drag coefficient. 

2. This analysis should be redone using data for a 
carbon-carbon nose and new data for ice erosion found in 
part 1 above. The mass loss for a carbon-carbon nose tip 
would tend to be far less than for the silaceous nose tip 
which is modeled in this thesis. Hence, mass loss due to 
general erosion would tend to have less effect than in this 
analysis . 
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It may be, however, that the mass loss due to ice erosion 
for the carbon-carbon nose is of the same order of magnitude 
as the mass loss due to ice erosion for the silaceous nose. 

If so, .the effects of ice erosion may be more important in 
the carbon-carbon nose than in the silaceous nose because 
the effects of mass loss due to ice erosion were masked in 
the silaceous nose by the substantial mass loss occurring 
throughout the reentry. 
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APPENDIX A 



POLYNOMIAL REPRESENTATION OF NOSE SHAPES 

The following pages contain graphs of the various nose 
shapes. Each graph contains a plot of the data points taken 
from the respective nose shape curves in Figure 10 with the 
aid of the HP-9874A digitizer and a plot of the polynomials 
which describe the data points. The polynomials also are 
displayed on the graph; the fifth-order polynomial applies 
for 0 £ x <_ 2.4 and the first-order polynomial applies for 
x > 2.4. 

In the analysis, the polynomial was used to describe the 
first 2.4 inches of the nose shape. The linear approximation 
was used to describe the nose shape from 2.4 inches aft of 
the nose tip to about 20 inches aft of the nose tip. The 
linear approximation to Figure A-l was used to describe the 
nose shape from about 20 inches aft of the nose tip (where 
the two linear portions intersect) to the end of the reentry 
body for all nose shapes. 

The profiles for nose shapes with ice erosion use the 
same equations to describe the shape. The coefficients are 
modified by substitution of a modified altitude into the 
equations for the coefficients displayed in Appendix B. 
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Figure A-l. Nose Shape at 250,000 Feet 
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Figure A-2. Nose Shape at 90,000 Feet 
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Figure A-3. Nose Sliape at 56,250 Feet 




in 

id 

i 

u 

z 



X 



(S3H0NI) A 



75 



figure A-4. Nose Shape at 40,100 Feet 
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Figure A-5. Nose Shape at 25,000 Feet 
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figure A-6„ Nose Shape at 15,500 Feet 



APPENDIX B 



POLYNOMIAL REPRESENTATION OF COEFFICIENTS 

The following pages contain graphs of the various nose 
shape coefficients defined in equations (49) through (52) . 
The first six graphs correspond to the coefficients used in 
equations (49) and (51) and apply for 0 <_ x <_ 2.4 inches 
where x is centerline distance from the nose tip. The last 
two graphs correspond to coefficients used in equations (50) 
and (52) and apply for x > 2.4 inches. 
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Figure B-l. Graph of Coefficient A(5) versus Altitude 
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ALTITUDE (FEET) 

Figure B-2. Graph of Coefficient A(4) versus Altitude 
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ALTITUDE: (FEET) 

Figure B-3. Graph of Coefficient A(3) versus Altitude 
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RLTITUDE (FEET) 

Figure B-4. Graph of Coefficient A(2) versus Altitude 
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ALTITUDE (FEET) 

Figure B-5. Graph of Coefficient A(l) versus Altitude 
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RLTITUDE (FEET) 

Figure B-6. Graph of Coefficient A(0) versus Altitude 
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ALTITUDE (FEET) 

Figure B-7. Graph of Coefficient B(l) versus Altitude 
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ALTITUDE (FEET) 

Figure B-8. Graph of Coefficient B(0) versus Altitude 



APPENDIX C 



POLYNOMIAL REPRESENTATION OF NOSE TIP RECESSION 

The following graph shows the evaluation of nose tip 
recession, the distance the nose tip ablates from the 
original location, as shown in Figure 10. As for the 
coefficients described in Appendix B, the nose tip recession 
was divided into two portions. A sixth-order polynomial 
describes the nose tip recession below 90,000 feet and a 
linear model describes the recession above 90,000 feet. 
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ALTITUDE (FEET) 

Figure C-l. Graph of Coefficient XZERO versua Altitude 



APPENDIX D 



COMPUTER PROGRAM EXAMPLE 

Appendix D is a sample of the computer program used to 
conduct the trajectory analysis. The computations were 
conducted on an HP-9845A computer which uses BASIC computer 
language. In general, the program is self-explanatory in 
that the names of variables are used whenever possible. 

Line numbers beginning with 1XXX describe the initial 
conditions. The counter K is zero until zero altitude is 
reached and then is equal to one . The counter L is zero 
during the first iteration and allows the establishment of 
initial volume, mass, and base area of the reentry body. 
After the first iteration of nose shape, L is equal to 1. 

The counter Z is zero until the altitude of the cloud is 
reached. The value of Z is changed to one during the 
iteration to find the mass loss due to the impact of the 
ice in the cloud. After the iteration is complete, the 
value of Z is two. If Z is set equal to two initially, the 
baseline trajectory with no ice ablation will be calculated. 

Line numbers beginning with 2XXX calculate the new 
position of the body as it descends. Line numbers beginning 
with 4 XXX calculate all of the major parameters needed for 
the calculation of pressure coefficients and the nose shape 
coefficients. Line numbers beginning with 5XXX calculate 
the actual nose shape, pressure coefficient, drag, volume. 
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etc. Line numbers beginning with 6XXX calculate the initial 
volume, mass and base area and values of drag coefficient 
and ballistic coefficient during the entire reentry. Line 
numbers beginning with 7XXX print the final values obtained 
during reentry and stop the program. Line numbers beginning 
with 8XXX calculate the mass ablated during passage through 
the ice cloud. Line numbers beginning with 9XXX calculate 
the final range and set the final altitude to zero so that 
calculation of other parameters at zero altitude is possible. 
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TRAJECTORY WITH ICE RELATION 
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7510 ! **************************^********************it*****’********i(’*'*** 

7511 ! THE FOLLOWING SECTION CALCULATES THE MASS LOSS DUE TO ICE ABLATION 

ALSO CHANGES THE NOSE SHAPE TO ACCOUNT FOR THIS MASS LOSS. 
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